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Abstract 

In the short wavelength Umit the BogoUubov quasiparticles of trapped Bose- 
Einstein condensates can be described as classical particles and antiparticles 
with dynamics in a mixed phase-space. For anisotropic parabolic traps we 
determine the location of the resonances and study the influence of the sharp- 
ness of the condensate surface on the appearance of chaos as the energy of the 
quasiparticles is lowered from values much larger than to values comparable 
with the chemical potential. 

I. INTRODUCTION 

The achievement of Bose-Einstein condensation of clouds of magnetically trapped al- 
kali atoms by evaporation cooling to temperatures in the 100 nano-Kelvin regime 
has revived the interest in the physics of weakly interacting Bose-condensed gases. As is 
well-known the bulk properties of such gases are rather well described by ideal gases of quasi- 
particles, first derived microscopically by Bogoliubov [Q, which are collisionless phonons at 
long wavelength and approach free particles at short wavelength. In the case of trapped 
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condensates of repelling atoms, to which we confine our discussion here, the quasiparticle 
description retains its usefulness. However, trapped condensates are spatially inhomoge- 
neous. They form around the minima of the trapping potential and can form rather sharp 
surfaces. The thickness of the surface layer, given by the healing length 0, can become 
very small compared to the radius of the condensate. Quasiparticles may therefore leave the 
condensate and reenter it after being refiected back by the trapping potential. Therefore 
their dynamics are more complex than in a homogeneous system. In Bogoliubov's approach 
single quasiparticles are described by wavefunctions, which are solutions of a set of linear 
wave equations. This is in complete analogy to the description of single particles in quan- 
tum mechanics by Schrodinger's equation. We know that it is extremely fruitful to examine 
the classical limit of the Schrodinger equation, which gives all of classical physics. In the 
same spirit we can examine the classical limit of the wave equations for the quasiparticles. 
In spatially homogeneous condensates this gives a simple theory of free quasiparticles with 
conserved momentum which are distinguished from free particles merely by their unusual 
relation between energy and momentum, and hence also between velocity and momentum. 

In the case of trapped inhomogeneous condensates recent work has shown |^ that the 
classical limit of the dynamics of quasiparticles becomes more interesting because they expe- 
rience forces both from the trap and from the condensate. In the case of isotropic traps, and 
hence also isotropic condensates, angular momentum conservation ensures the integrability 
of the quasiparticle dynamics. This can be used to construct WKB solutions of Bogoliubov's 
wave equations for this case 0. In the experimentally more relevant case of axially sym- 
metric traps integrability of the quasiparticle dynamics is lost and numerical studies 
indeed have shown a generally mixed phase-space in this case. Detailed analytical work was 
performed for quasiparticle energies E much smaller than the chemical potential fi. Here 
two integrable limits were identified: (i) the phonon limit, where the phonon-like quasipar- 
ticle is confined to the interior of the condensate and is specularly refiected back when it 
strikes the surface, and (ii) a surface-particle limit, where the motion of the single-atom-like 
quasiparticle consists of rapid small-amplitude oscillations between the repelling main bulk 
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of the condensate and the potential wall of the trap and a slow secular motion along the sur- 
face of the condensate. The numerical examination of the quasiparticle dynamics at larger 
energies (and for the trap anisotropy of the experiment [|1|) revealed ||^ a strong chaotic 
component at E = fi and a curious ' quasi- integrable' regime at E ^ fi, where the dynamics 
in phase-space clusters around the tori corresponding to single atom motion in the trapping 
potential, however with small-scale chaos superimposed. 

In the present paper we wish to study this large energy regime more closely, both ana- 
lytically and numerically. In the spirit of Chirikov's pioneering work on the onset of chaos 
(reviewed in 0) we ask 'where are the resonances' and answer this question by developing 
the first steps of the classical perturbation theory for the quasiparticle dynamics at large 
energy. We also examine the influence of the thickness of the surface on the appearance 
of chaos. In the Thomas-Fermi approximation |Q this thickness is neglected, leading to a 
spatially discontinuous effective force on the quasiparticle. This violates assumptions of the 
KAM theorem and turns out to be the reason that noticeable small-scale chaos survives 
even at very large energies, where the system should be close to the integrable limit of in- 
dependent atoms in the trap. Simulations with different boundary layers substantuate this 
hypothesis. For a condensate with boundary layer we study the transition to chaos as the 
energy is lowered to values comparable with the chemical potential. 



II. EQUATIONS OF MOTION 

In the present section we give a brief derivation of the relevant equations of motion. The 



starting point is the Gross-Pitaevskii equation [jTO| 

thij{x,t) = |-|;^V2 + Uix) + Vomx,t)\^'^ ^{x,t) (2.1) 

describing the macroscopic wavefunction ip oi a Bose condensate at temperature T = in 
a trap with potential U and interaction Vq = Anh^a/m, where a is the s-wave scattering 
length. We shall consider only the case of repulsive 2-particle interactions a > 0. The trap 
potential is taken as axially symmetric and parabolic 



U{x) = '^{x' + y' + \z'). (2.2) 

The ratio between axial and radial trap frequency is \/~X. The particle number = 
/ d^x\il){x,t)\'^ is conserved and assumed fixed. The equilibrium state of the condensate 
with chemical potential is a solution of ( p.lD with il){x^t) = e~^^^/^il)Q{x) which satisfies 

liMx) = \ + U{x) + V,\M^, t) |2 1 ^,{x) (2.3) 

with vanishing axial angular momentum 

= - / d^xr,{x){x X V)^^{x) . (2.4) 

The small oscillations around the equilibrium state are waves of the form 

i,{x, t) = e-*^*/^ (Mx) + Vi^, t)) (2.5) 



with {p considered as small and satisfy a linearized version of eq. (|2.1| ) 

ih^{x,t) = l^-^V^ + U{x)-fi + 2Vo\Mx)\^'^vix,t) + Voij^o{x)^*ix,t). (2.6) 

Usually the small oscillations are considered in second quantization (where (|2.6| ) becomes 
an operator equation, but retains its form), which brings out their particle-like properties, 
and they are then called quasiparticles. Here we shall bring out the particle properties in a 
different way, by considering the 'ray-optics' hmit of eq. ( |2.6| ) [0. Formally, this is done by 
means of the ansatz 

ifix, t) = (ao(a;, t) + 0{h)) e*^(-'*)/^ - %{x, t) + 0{h)) e-^^(-'*)/^ (2.7) 



and considering the terms arising from eq. (|2.6| ) formally order by order in %. The sym- 
metrical form of the ansatz (|2.7| ) allows us to restrict the sign of dS/dt by choosing 
dS/dt = —E < 0. flo and Bq can be interpreted as semiclassical amplitudes of the quasipar- 
ticles and their antiparticles, repectively. Writing p = VS" for their momentum one finds to 
zeroth order coupled algebraic equations for the amplitudes oq, &o 



^ ^ ao{x,t)^ 



^ -Voipfix) en Fix, p)- 



dS{x,t) 
dt J 



V 



(2.8) 



where 



eHF = ^ + U{x) - /i + 2Vo\Mx) 
2m 



(2.9) 



is called the classical Hartree-Fock energy. 

Eq. ( |2.8D implies as solvability condition a quadratic equation for dS/dt which reduces 
to the Hamilton- Jacobi equation 



^ + H{x,VS) = 



(2.10) 



with 



H{x,p) = Jelp{x,p) -VilMxW 



(2.11) 



where by the sign convention on dS/dt the positive branch of the square-root must be 
taken. Eqs. ( p.lO| ), (|2.11|) now describe the classical dynamics of the quasiparticles comple- 
mentary to the waves with dynamics ( p.6| ). The conservation of energy H and axial angular 
momentum of the quasiparticles can be taken care off by separating 



S{x,t) = So{p,z)-Et-L, 



(2.12) 



where p, 0, z are standard cylinder coordinates. Eq. ( |2.10| ) then reduces to 

^ dSo dSo^ 



H \p,z, 



E 



(2.13) 



dp dz ^ 

with \ijjoix)\'^ = \Mp,z)\'^ andp2 = p'^^+p^^ Ll/ p^ in eqs. (|2j), (|2Tl|) . Eq. (U) fixes only 
the ratio of the amplitudes 6o and oq, which, for fixed E and L^, becomes 

\E^ + V^\Ux)\'f"-E^"' 



(2.14) 



_{E^ + vimx)\^y'^ + E_ 

In order to determine the absolute values (apart from an arbitrary space-independent factor) 
one has to go to the next order in the expansion where one finds as a solvability condition 
the conservation law 



^(|aor - l&oH + ■ (|aoP + \bo\')VS = 



(2.15) 



In the following we shall be concerned only with the classical quasiparticle dynamics de- 
scribed by eq. (KW) with the Hamiltonian (|2Tl|), 



III. PERTURBATION THEORY FOR LARGE ENERGY 

The single-particle interaction energy with the condensate is of the order of the chemical 
potential fi and is only a small perturbation for energies E ^ /i. Expanding to second order 
the Hamiltonian takes the form 

H = eo{x,p) + 2Vo\M^)\' - ^ (3.1) 

with 
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eo = 1^ + j^oi^' + y' + A2;2) - ^^ . (3.2) 

In the following we shall restrict our discussion mainly to the case of vanishing axial angular 
momentum = 0, but at the end we also present some results for 7^ in order to 
assess to what extent the case = already captures the typical behaviour of the system. 
Experimentally modes with = or 7^ can be excited depending on the symmetry of 
the excitation mechanism. 

For Lz = the dynamics is restricted to a plane containing the 2;-axis, which can be 
taken as the (x, 2;)-plane y = 0, Py = 0. The action angle variables of the unperturbed 
harmonic motion in the trap are 



X = sinOr, , = J 2muJo\x cos 9^ 

V mujQ ^ 



(3.3) 



— — sin 9z , Pz = \l 2m\/~XuJo\z cos 9z 
ray Acuq 



with 
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eo = ujqIx + \^uJolz - /i . (3.4) 

To express the perturbed Hamiltonian in action-angle variables we need the Fourier coeffi- 
cients of the condensate density and its square 



(2vr)^ 



JO 



(3.5) 



p2,^(i.,i,) = jj— r rde,e-'^''-+-'^^\M^,o,z)\^ 



(27r)2 Jo 

The canonical transformation {6, 1) —>■ (0, J) with 



^0 7;;^£ + vAn 



2Vo 1 dpi^n{Jx, Jz) ^i(W^+ne,) 



(3.6) 



'^0 ^ i{i+y/Xn) dJx 
and analogous for I^, 0^ removes the angle-dependence in the first order perturbation term 
and we are left with the second-order Hamiltonian 

H = uoiJx + VXJ,) - /i + 2VoPooiJx, Jz) + — J2 ^An(4, J,)e*(^<^^+"<^^) (3.7) 

where 

= 4 y ^P£-p,n-r- Pp,r .g g^ 

Eq. ( p.7| ) differs from ( p.l| ) only by terms of higher than second order. For irrational there 
are no resonances to zeroth order in the interaction, but in first order isolated resonances 
^x{Jx, Jz)P + ^z{.Jx, Jz)R = appear with integers P, R and the perturbed frequencies 

<9poo 



Vtx = cuo + 2Vo 



dJx 

(3.9) 

dJz 

Their distances in frequency space have to be compared with their widths, given by 



where 



Yl \Kpn{J.,J.)\ 



v^2 p 8'^ POO p 

^'''=^ '"■^-"■'^ P. = P, P^ = R. (3.11) 



Mpn{Jx,Jz) P^ + PI 

If resonances overlap somewhere in phase space Chirikov's criterion tells us that we should 
expect chaos there. Unfortunately, eqs. (|3^), ( |3.8| ) are too difficult to evaluate analytically 
for a realistic equilibrium solution ^0(2^) of the time- independent Gross-Pitaevskii equation. 
However, these equations are still useful to understand some qualitative features of numerical 
simulations of the complete Hamiltonian dynamics at large energies. 

Let us discuss in particular the case of large condensates to which the Thomas-Fermi 
approximation applies, in which the spatial derivative terms in eq. ( p.3|) are neglected. The 
condensate density is then given by 



Vq I 

where ^(x) is the Heaviside step-function. As |?/'o(a;)p is an even function of x, y, z the per- 
turbation amplitudes are nonvanishing only for even £, n. Therefore Kpp^ in ( |3.10D must be 
replaced by K2P2r- The Thomas-Fermi approximation of the condensate density (|3.12|) has 



a discontinuous first-order derivative at the surface. This means that the Fourier coefficients 
pe^n and p^^^ at large \n\ fall off like 1^*1 |^|~^ and \n\~^, respectively giving the 



(P, -R)- resonances for large |P|, |P| widths, which according to eq. (|3.8| ), fall off only like 
|P|~^, |P|~^. This estimate results from the second term in eq. ( p.8|) , which is predicted 
to fall off only like |^|^^. On the other hand the number of large-order resonances 

Qr^/Qz = —R/P scales like |P|, |P|. Therefore, barring non-generic cases where the l-ft'piji 
are small for some exceptional reason, one expects large-order resonances to always overlap 
in the Thomas-Fermi approximation, i.e. the tori of the free harmonic oscillations in the 
trap will typically all be broken. By the same arguments a condensate density with M 
smooth derivatives and a discontinuous M + 1 order derivative will give rise to resonance 
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widths scaling like \P\~^~^, \R\~^~^ and tori with large |P|, \R\ can exist if M > Mc = 2. 
The critical order of smoothness for general Hamiltonian systems with / degrees of freedom 
was determined by Chirikov [^] as Mc = 2/ — 2. 

In fig. 1 a Poincare surface of section of the dynamics a.t E = 20/i, Lz = 0, taken at 
X = and plotted in the (6'^, l2)-plane is shown for the experimentally realized [HI value 
A = 8. The action variable I^ is plotted in units of With the exception of tori at 

small lz and at large I^ ~ {E + fi)/ ^/Xluq all tori in fig. 1 are broken. The survival of the 
exceptional unbroken tori can be understood from the fact that for ^ or —>■ (with 
Iz ^ {E + fj,)/^/XuJo) all the coefficients K^^n approach with the exception of K^^q or Ko^n, 
repectively, which cannot influence appreciably tori with frequency ratios n^/^lz — 8^^^^. 

In fig. 1 a number of resonances with frequency ratios in the neighborhood of Q^/^z = 
8~^/^ can be discerned. Their frequency ratios are given on the right hand side of the 
graph. To understand the position where these resonances occur we consider poo in the 
Thomas- Fermi approximation 
4 



Poo 



TcWo Jo 



J d(f)x j d(j)z (^/i — ojqJx sin^ 0^. — ^/XujqJz sin^ Cf 

(fi — uqJx sin^ <^a; — VXujqJz sin^ (pz) (3.13) 



and evaluate the first-order frequency shifts 



Au,,z = (3.14) 



from the integrals 

I'Tt/2 riv/2 



AWa; = 5- / / #2 siu^ 4>xQ{\^ " ^oJ x sin^ - V Xu^qJ z sin^ 0z) 

TT^ JO JO 



(3.15) 



Aa;^ = — / ci0a; / d<\)z sin^ (j)zO{fi - uJqJx sin^ 0i - v XuqJz sin^ 0^) . 

Tl^ Jo Jo 

It is manifest that AuJx,z are negative. Near the bottom of fig. 1 the action Jz is small, while 
Jx — {E + fi) / ujo — \fXJz is large. The Heaviside function in eqs. ( |3.15| ) therefore restricts 
sin^ (px to small values of the order of fi/uoJx = 0{^/ E) while -\/A sin^ 0^ is not so restricted. 



Therefore lAuJ < |Au;^| near the bottom of fig. 1 and ^ 7^ ^ - ^ > holds 

there. In the upper parts of fig. 1 Jx is small while Jz — {E + j^)/ ^/Xujq — J^/ \/~\- The 
situation is therefore reversed and we obtain ^^^^ — ^ < by the same argument. 

The periodic orbit z = 0, = which has = forms the lower border of the range 
of Jz- For this case the frequency shifts Awa^, Aa;^ are easily evaluated from eq. (|3.15| ) with 
the result 



= -O (oJoi^i/uJoJxf^) 



(3.16) 



2 XuJoli 



TT V Jx 

indicating that the fixed point is stable for E ^ fi. The arguments of the previous paragraph 
indicate, furthermore, that \AuJx\ is minimal for this case, while \AuJz\ is maximal, leading 



to a maximal value of the ratio 



(3.17) 



-^2 /max 7l\l XtOoJx \ 71 \j E + fl J 

where we used + yU ~ tuo + VX^^oJz with = in the last estimate. For E/ fi = 20 this 
gives (Qx/^z) — 0.402, which is just barely larger than the ratio 0.4 of the 2:5 resonance 
visible near the lower border of fig. 1. Similar arguments apply to the upper border of the 
range of Jz which is formed by the periodic orbit x = 0, px = with vanishing Jx for which 



Aoox = -IJ^^ > = -V^O {{f^/VXio.Jz)'/') . (3.18) 



- ' • (3.19) 



Thus 



^^z/min V 7^]] E + 

For E/fi = 20 this yields (^^a;/^z)min — 0.304 which is smaller than the ratio 0.333 of the 
1:3 resonance visible near the upper border of fig. 1. 

In general we can conclude that for fixed energy E we should expect to see the strongest 
resonances in the interval 
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which are, for A = 8, E/ ^ = 20, the resonances 1:3, 2:5, 3:8, 4:11, 5:14, 6:17 and are indeed 
all seen in fig. 1 with the exception of the 5:14 resonance. In its place a chaotic region is seen. 
The resonance values of the actions following from the perturbation theory are indicated on 
the right-hand side of fig.l. The 4:11 resonance predicted close to the missing 5:14 resonance 
is nearly completely destroyed and barely separated from the chaotic region replacing the 
5:14 resonance. 

Let us now consider the effect of the sharpness of the boundary layer of the condensate. 
Its thickness is given by the healing length 

It enters as an independent parameter whose ratio to the radial Thomas- Fermi radius Ptf = 



2fi/muQ is determined by the value of the chemical potential, i.e. by the number of particles 
in the trap, in/ Ptf = fiujQ/2fi. As example we consider in/ Ptf = 0.1. The condensate with 
boundary layer is modelled by joining 

\M^)\' = ^~\'?'"'"' 0{pTF-iH-r) 

+ exp(-ao - ai(r - ptf + ^h) - a2{r - Ptf + ^nf (3.22) 
-a^{r - PTF + ^Hf )0{r - Ptf + (-h) 

continuously with continuous three derivatives at r = p^F — ^h, thereby fixing oq, . . . , 03. 
Here the abreviation r = y/ p'^ + Xz'^ is used. The degree of smoothness thereby introduced 
should be sufficient to see smooth KAM tori. The model condensate (|3.22|) is of course not 
a solution of the Gross-Pitaevskii equation, but does not differ qualitatively from solutions 
taking into account the boundary layer [|1^]. For our present purposes it is therefore per- 
fectly acceptable while avoiding unnecessary complications. In fig. 2 we compare Poincare 
surface of sections in the same format as in fig.l for the Thomas- Fermi approximation and 
the smooth model condensate with Ih/ptf = 0.1 at energy E/p = 100 and A = 8. It 
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can be seen that with the smooth boundary layer included most of the broken tori of the 
Thomas-Fermi condensate seen in the upper part of fig. 2 become smooth, as seen in the 
lower part of fig. 2, however with ripples on the tori still indicating the presence of the nar- 
row boundary layer. If the boundary layer is narrowed further (not shown) these ripples 
become stronger and develop sharp cusps. The chaotic band at small actions survives 
even for very large quasiparticle energies E/ fi. Resonances within this band have a com- 
paratively large width for two reasons. For one the factors 1/MpR are large, because Iz is 
not far above the value Iz = n/uooV^ where Mp^ peaks (with a logarithmic singularity in 
Thomas- Fermi approximation). For larger values of Iz, of the order of E /ujQ\f\, like for the 
6:17 resonance or values oi P : R close to y/X, l/MpR is much smaller by a factor of the order 
of {jji/EY. Second, the interaction coefficient \Kpji\ has a resonantly enhanced contribu- 
tion —4:{dpoo/dJz)p2P2R{P/R+ \/A)~^. The last factor in this expression favors resonances 
P : it! close to the ratio VX. The existence of the narrow boundary layer of the conden- 
sate furthermore leads to appreciable interaction coefficients \Kpii\ even for comparatively 
high-order resonances whose overlap may be responsible for the formation of the chaotic 
band. The requirement that both factors Mp^ and \K2P2r\ must be large may explain the 
appearance of a band of actions Iz with appreciable resonance overlap somewhat above the 
value Iz = p/ljoVX- 

Dynamically the chaotic band is related to orbits in a layer of actions Iz ^ ii/y/X/ujQ 
which are just sufficient for the quasiparticle to hit or miss the condensate at random as 
it oscillates back and forth in x-direction. This mechanism is able to introduce sensitive 
dependence on initial conditions and exists, if only in a narrow band, even at very large 
energies. For just shghtly smaller actions Iz the quasiparticle has to pass the condensate 
twice in each period of 6^ and the tori are smooth even in Thomas-Fermi approximation, 
as can be seen in fig. 1 for Ef/i — 20. A similar mechanism for a chaotic band should 
actually exist also for the oscillations in the z-direction, the short axis of the condensate- 
ellipsoid. However the instabilty in this case for action I^ ^ n/\/\/uJo seems to be much 
less pronounced (which is plausible, because the perturbation by the condensate should be 
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weaker along the short axis) and is not visible in the numerical data. 

Finally we consider numerically the transition to chaos in the condensate with boundary 
layer Ih/ptf = 0.1 as the quasiparticle energy is lowered. In fig. 3 the same Poincare surface 
of section is plotted as in figs. 1,2 (and again for A = 8) but in the (2;,p2)-plane rather than 
the (^2, /2)-plane. The two plots in the upper row are for E/ ji— 100 and E/ ji — 20, those in 
the lower row are for E = \Q and E = 2 respectively. The plots for E = 100 and E = 20 can 
be compared with fig.2 and fig.l respectively, but fig.l is, of course for vanishing Ih only. 
The shghtly rounded Thomas- Fermi surface at z/ptf = A~^/^ is visible in these graphs. 
Based on these and similar plots the transition to chaos can now be roughly described as 
follows: The chaotic band in a layer of small actions ^ p/VXcuo existing even at very 
large energies (see the plot for E/p — 100) becomes wider as the energy is lowered (see 
E/p = 20). Then, at a critical energy E/p = 13.9 which we discuss further below, the 
periodic orbit along the long axis of the condensate-ellipsoid becomes unstable and a second 
chaotic region in phase-space surrounding the unstable orbit z — — is created. The 
inner and outer chaotic regions then rapidly grow together as the energy is lowered further 
(see E/ jjL — 10) and finally fill up most of the accessible regions of phase-space (see the last 
plot with E/ p — 2). 

Apart from a continuous widening of the chaotic regions the main event in this transition 
to chaos is the appearance of the instabilty of the periodic orbit along the long axis of the 
condensate-ellipsoid, which for A > is the x-axis. In Thomas-Fermi approximation this 
orbit has the time-dependence for \x\ < ptf 




(3.23) 



Px = iT^^oPTFy E / p — 1 cosh(a;o^ ~ <^°) 




and for \x\ > ptf 
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(3.24) 

Px = muJoPTFy E/ fi + 1 cosh(u;ot — f°) 



with ip° and (p'° suitably adjusted by continuity at |x| = ptf- To determine the energy where 
this orbit looses its stability we can perform a linear stability analysis for small perturbations 
z, pz away from this orbit, which satisfy 

mz = pz, Pz = \ {I — 29{\x{t)\ — Ptf)) ^^qZ . (3.25) 

The growth of a perturbation of the periodic orbit is determined by the monodromy matrix 
M for a half-period T/2 



\ ( z^^)/p^- ^ 



M 



TF 



\^pfymUQPTF j 



(3.26) 



z{TI2)Iptf 
yPz{T/2)/muJopTF j 

The Hamiltonian form of the dynamics ensures that Det iVf = 1. The stability condition 
for the perturbations in z-direction then becomes 

\TrM\<2. (3.27) 

With a little algebra it is straight-forward to evaluate M and its trace thereby reducing 
O to 

I cos [yXijJi^t'^ cosh (^-v/Acjo^i) I ^ 1 (3.28) 

Here t\ and ^2 are the total lengths of the time-intervals during each half-period where 
|x| < Ptf and |x| > ptf, respectively. They are given by 

ti = — Artanh S^, ^2 = — Arctan W — . (3.29) 
a;o V E a;o \ p. 

For A = 8 stability is lost, according to the criterion (|3.28| ) for Ejp = 14.4 which is in rea- 
sonable agreement with the already quoted value Ejp = 13.9 determined from the numerical 
simulation of the complete quasiparticle dynamics. 
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IV. CONCLUSION 



In the present paper we have studied the classical limit of the dynamics of the quasipar- 
ticles in a spatially inhomogeneous harmonically trapped weakly interacting Bose gas. As 
is well known many properties of a Bose gas are determined by its quasiparticles. For suffi- 
ciently large energies E /Tiujq ^ 1 a classical description of quasiparticles should be possible. 
We have accordingly concentrated our attention on the large energy regime. For quasipar- 
ticle energies large compared to the chemical potential a perturbative expansion becomes 
possible. The dynamics then looks like that of an anisotropic harmonically bound particle 
which is perturbed by the presence of a weakly repelling condensate localized around the 
center of the trap with a rather sharp surface. The appearance of an infinitely sharp surface 
in the Thomas-Fermi approximation leads to the break-up of all tori with the exception 
of those in the immediate phase-space neighborhood of the periodic orbits along the main 
axes of the ellipsoidal condensate. Therefore in this case chaos exists in phase-space for 
arbitrarily large energies E / ^. Still resonances can be identified even in this case and their 
actions are well described by perturbation theory. We have compared this non-smooth case 
to that of a condensate with a smooth (up to third derivatives) but narrow boundary layer. 
In real condensates the thickness of the boundary layer is determined by the two-particle in- 
teraction and the number of particles |^ . In that case smooth KAM tori exist at sufficiently 
large energies, but they are rippled by the infiuence of the boundary layer. Furthermore it 
turns out that an appreciable region of chaos in phase-space persists even to large energies. 
It is related to orbits along the long axis of the condensate-ellipsoid which are sufficiently 
perturbed in the direction of the short axis to sometimes miss and sometimes hit the con- 
densate in a random way with sensitive dependence on small perturbations. Finally we have 
studied how large scale chaos appears as the energy is gradually lowered to values of the 
order of the chemical potential. The instability of the periodic orbit along the long axis 
of the ellipsoid was found to play a major role in this transition. It is connected with the 
appearance of a second inner chaotic region at lower energies which joins up with the chaotic 
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band existing also at large energies. 

Our discussion so far has been restricted to the special case = 0, where the classical 
quasiparticle dynamics is confined to a plane in configuration space containing the z-axis. 
One may well ask to what extent this motion already captures the typical behavior of the 
system. To examine this question at least numerically we have generated Poincare surface 
of sections for Lz ^ both for fixed E and varying Lz or for fixed Lz and varying E. The 
surfaces of section are always taken at the value p = po where the effective potential in 
radial direction, which includes the centrifugal barrier, has its minimum. In fig. 4 we present 
a series of Poincare surfaces of sections for Lz fixed at a rather high value (in units of h) 
Lz = p/^o and the same values of E and also the same thickness of the boundary layer as 
in fig. 3. 

It is apparent that taking 7^ the mirror symmetry with respect to the z-axis and the 
j9^-axis which is present at = is lost and replaced by a point symmetry with respect to 
the origin 2; = = 0, which is, of course, to be expected. However, apart from this obvious 
difference the qualitative behavior displayed in fig.4 is the same as in fig.3. As was shown 
in the energy and the angular momentum must satisfy the inequalities 

E + p>uoLz if E>p (4.1) 
E > {uoLzY/Ap if E < p. 

For E > p and \Lz\ > 2p/uJo there exists a domain {ojqLzY j-ip > E > uqLz — p where 
the motion of the quasiparticles occurs outside the condensate in the Thomas- Fermi limit, 
and is therefore integrable. Thus an additional transition from nonintegrable to integrable 
dynamics occurs for very high angular momentum \Lz\ > 2p/ujQ as the energy is lowered 
from above to below E = {ujqLzY / 'ip. 

Let us finally turn to the quantum mechanical implications of the classical dynamics 
we have studied. Such implications exist both for the wavefunctions and the energies of 
single quasiparticle states which are of course closely connected. As the classical analysis 
is done in phase space it is most convenient for a discussion to use the Husimi distribu- 
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tion Q{p,Pp, z,Pz) — I < alip > p where jo; > is a coherent state of the harmonic oscillators 
defined by the free trap. For each quantum state \ip > the function Q is a positive quasiprob- 
ability on the phase space, which due to the overcompleteness of the \a >, contains the full 
information on >. 

Let us now turn to the Poincare surface of sections displayed in fig. 3 and discuss the 
corresponding quantum states via their Q- functions: In fig. 3a the Q- functions at p = 
will be spread out along the tori with (in the limit Tiujq/ ^ — > 0) narrow peaks on the tori. 
The number of such states is semiclassically given by the phase space volume of the energy 
shell of thickness in units (27r^)^. The chaotic layer visible in fig. 3a may correspond 
to several quantum states (depending on the contribution of this layer to the phase space 
volume of the energy shell in units (27r^)^) which are all spread out along the layer but with 
wave functions oscillating wildly in the corresponding configuration space so as to satisfy the 
orthogonality condition. If indeed several such states exist in an energy shell of thickness 
/S.E their energies will tend to repell each other. The energies corresponding to Q-functions 
centered on tori will show no such repulsion. 

Essentially the same discussion applies to fig. 3b, 3c where the chaotic layer occupies a 
larger fraction of the phase-space volume and finally merges into a single domain. In the case 
of fig. 3d, where essentially a single chaotic domain survives two different types of quantum 
states are possible: The first possibility is that the quantum states are all spread out over the 
chaotic domain with strong oscillations of their wave functions in the corresponding domain 
of configuration space to satisfy orthogonality and accompanying strong level repulsion. The 
second possibility is that the quantum states show dynamical localization with respect to 
the action variable 7^. In this case the quantum states still localize around different values of 
Iz due to a coherent interference effect akin to Anderson localization. In this case the wave 
functions can be orthogonal without energy level repulsion. Which of these two possibilities 
is realized depends on the strength of the chaotic change of the action variable 1^. If this 
change is sufficiently large, then the first possibility (extended states with level repulsion) 
will be realized. If however Iz changes sufficiently weakly in a diffusive way such that its 
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variance satisfies < A/f >= Dt with a diffusion constant D, then the second possibihty may 
be reahzed. Then a fundamental estimate [|1^] can be made of the locahzation length of 
Iz according to S,1 ~ Dhp{E) where p{E) is the density of states at energy E. Dynamical 
localization then is predicted to occur as soon as C,z <^ Wz where Wz gives the width of the 
chaotic domain in J^. Thus D ^ /hp{E) is required. At present no reliable estimate of 
D exists, either analytically or numerically, so the question which of the two behaviors will 
occur in a given trap at energies E p must be left open here. Numerical computations of 
quasiparticle energies have been performed for certain trap parameters and avoided crossings 
have been seen in the results for the regime ~ /i (see e.g. [0)' t>ut a systematic study of 
the level statistics has not yet been performed for this case. 

However, the level repulsion typical for classically chaotic systems has already been built 
into a quantum theory of damping of the low-energy collective modes due to scattering of 
thermally excited quasiparticles at energies of the order of the chemical potential . This 
should only be a beginning. If the precedence of mesoscopic systems is any hint, we may 
expect that more properties of trapped Bose condensates will turn out to be infiltrated by 
chaos in the future. 
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FIGURE CAPTIONS 

Fig. 1. Poincare surface of section of the quasiparticle dynamics. The cut is taken at con- 
stant energy at x = 0, > for energy E/fi = 20 and A = 8. The cross-section is presented 
in the (6*2, l2)-plane, where y/Xuolz = p1/2m + {m\uj1/2)z^ , arctan6'2 = m\/~XujQz/pz- The 
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action variable is given in units of 2/i/a;o. Frequency ratios ft^ '■ of resonances are given 
on the right hand margin together with the resonance actions determined from perturbation 
theory. For resonances existing in doublets related by symmetries only one member of the 
doublet is shown for clarity. The 5:14 resonance could not be detected numerically. 
Fig. 2. Poincare surface of sections as in fig.l for E/fx = 100 for the dynamics with 
the condensate in Thomas-Fermi approximation (upper part) and for the condensate with 
boundary layer of thickness Ih/ptf — 0.1, continuous up to and including the third derivative 
(lower part) 

Fig. 3. Poincare surface of sections as in figs. 1,2 but plotted in the {z,Pz)-plsine. z,Pz are 
plotted in units '^2fj,/mwQ, ^/2mJI. The value of E/fj, is 100 (upper left), 20 (upper right), 
10 (lower left), and 2 (lower right). 

Fig. 4. Poincare surface of sections as in fig.3, but for finite angular momentum = /j^/ujq. 
Again the values ol E/ji are 100 (upper left), 20 (upper right), 10 (lower left), and 2 (lower 
right) . 
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